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Abstract 



The self-consistent spatial distribution of particles of Galactic dark matter 

is derived including their own gravitational potential, as also of that of the 

visible matter of the Galaxy. In order to reproduce the observed rotation 

curve of the Galaxy the value of the dispersion velocity of the dark matter 

particles, (w 2 )*/ 2 , should be ~ 600 km s -1 or larger. 
PACS numbers: 95.35 +d, 98.35 -a, 98.35 Gi, 98.62 Gq, 98.35 Df. 
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More than 20 years ago, it was suggested |]] that weakly interacting particles of non-zero 
rest mass which decouple from radiation and matter early after the Big Bang would form 
an invisible gravitating background of dark matter (DM) around galactic systems. Even 
though at that time the only available candidates for these particles were the neutrinos 
of the muon and electron flavors, the idea itself became the paradigm under which the 
newly discovered particles like the tau-neutrino and newly hypothesised particles within 
the context of possible physics beyond the Standard Model of particle physics could be 
incorporated. Also, during the latter half of the intervening twenty years, we have witnessed 
a tremendous growth in the experimental effort towards direct detection of these particles 
in the laboratory. The experiments are aimed at observing the effects of the impact of 
mainly the more massive candidate particles of DM with targets maintained at cryogenic 
temperatures which facilitate the observation of the tiny amount of energy deposited in the 
process against the background generated by internal and external radioactivity and by the 
cosmic rays. These developments are reviewed in detail by Trimble [Q] , Primack, Seckel and 
Sadoulet 0, Caldwell @, and Price § . 

The interpretation of these experiments to derive constraints on the properties of the 
unknown particles constituting a halo of dark matter in and around the Galaxy requires 
assumptions about the density and spectrum of velocities of the DM particles in the solar 
neighbourhood. These parameters have been obtained thus far by describing the DM halo 
as a single component isothermal sphere which is truncated at a particular radius ||. The 
normalization for the density of DM particles comes from an analysis originally suggested 
by Oort |7j in which the observed spatial- and velocity distribution of stars near the solar 
system indicate a DM density of ~ 0.3GeVcm -3 in the solar neighbourhood; Bahcall || 
gives a detailed account of this procedure. The 3-dimensional dispersion velocity of the DM 
particles, (w 2 )™, has not been determined, however. It is customary to take recourse to the 

virial result pertaining to an isotropic isothermal sphere [|| and set (v 2 )^ = \f\ Ooo ; where 
Gqc is the asymptotic value of the circular rotation speed. Since 6^ for the Galaxy is not 
known, the usual practice is to assume that the rotation curve of the Galaxy [|T^|n],[T^] , Q(R), 
is flat from R ~ 5 kpc out to R 3> -Ro ~ 8-5 kpc (here and below R denotes the galactocentric 
distance in the plane of the Galaxy, Rq being the sun's position), and set Goo ~ ©(-Ro) ~ 
220 km s _1 , the rotation speed near the solar system. This yields (v 2 )^ ~ 270 km s _1 , which 
is the value usually assumed in most studies of issues related to Galactic DM. However, as 



noted in the recent review by Fich and Tremaine ||12|| , "Much of the data indicates that the 
rotation curve continues to rise beyond Ro" . Thus the estimate (t> 2 ) 1//2 ~ 270kms _1 derived 
by assuming Boo = @(-Ro) is uncertain. Moreover, the assumption of a pure isothermal 
sphere for the description of the dark matter halo neglects the possible deviation from 
spherical symmetry induced by the disk-like distribution of the visible matter. 

Keeping these points in mind, we focus attention on the observed rotation curve of 
the Galaxy, and develop a theoretical framework, the salient features of which are: (a) A 
model for the Galaxy comprising of visible matter and particles of DM with a self- consistent 
inclusion of their gravitational interactions, and (b) Departure from spherical symmetry due 
to the disk-like distribution of the visible matter which will be treated as axially symmetric. 
The quantity (v 2 ) 1 ^ 2 appears as a free parameter in our framework and is determined by 
comparing the theoretical rotation curve with the observed data. 

We adopt well-established models to describe the density distribution of the normal 
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visible matter and the resulting gravitational potential. In this Letter we present our results 
for a two- component model of the visible matter consisting of a spheroidal bulge ]9 ,13,14 
with density p s (r), and an axisymmetric disk |TJ| with density pd{R, z): 



Ps( T ) = ^73/2 W 



p d (R,z) = '^-e-( R - Ro V R *e-W h (2) 

1 /2 

where r = (i? 2 + z 2 ) , and S = /f^ Pd(Ro, z )dz is the disk surface density at the solar po- 
sition, z being the vertical distance from the plane of the disk. The values of the parameters 
are given by |i"3| , |14| a = 0.103 kpc, Rd = 3.5 kpc, h = 0.3 kpc, andp s (i?o) = 7xlO _4 M pc~ 



(Note that the rotation curve in the outer regions of the Galaxy is relatively insensitive to 
the spheroid parameters). There are conflicting reports on the value of S : Whereas Kuijken 
and Gilmore |L4] (KG) suggest S ~ 40M Q pc -2 on the basis of data on ~ 512 K-dwarf 



stars, Bahcall et al ]r5| ) |TE| in their reanalysis of essentially the same data suggest a number 
for Eo which is about twice as large. In our calculations we consider values of So in the range 
(40-80) M pc -2 . The estimate of the local surface density of the Galactic disk due to the 
identified matter such as visible stars is ~ 48 ± 8 M Q pc~ 2 . Thus Bahcall et al's kinematical 
estimate of So seems to indicate the presence of a substantial amount of unseen matter in 
the Galactic disk, whereas KG's estimate is consistent with no disk dark matter. (Note 



that analyses of Refs. |14" , |l5 ,16| are all based on 1-dimensional solutions to the Boltzman 
equation, which, in the given situation, are strictly valid for an infinite disk only). In any 
case, the dark matter associated with the disk is likely to be dissipational in contrast to 
that constituting the extended halo which would be collisionless and non-dissipative. We 
are concerned with this latter type of dark matter in this paper. We use the conventional 
nomenclature "visible" to describe effectively the total matter associated with the disk and 
write the total visible matter density, p v , as p v = p s + pa, the corresponding potential being 

= $ s + $d. The expressions for the potentials $ s and $d corresponding to the chosen 
forms of p s and pd are given in Refs. [P|, [T3] , |T4]| . 

Now, for the DM component, the exercise is to calculate the distribution of the DM 
particles by self-consistently including the effects of the self-gravitation of the DM particles 
themselves and the potential due to the total visible component specified above. The proce- 
dure we follow is analogous to the one developed earlier |T7| with this difference that we now 



have to contend with the axial symmetry of the potentials. Since the DM particles obey the 
steady-state collisionless Boltzmann equation, the assumption of Maxwellian phase-space 
density allows us to write the spatial density, p DM (R, z), of DM as 



p DM (R,z) = p DM (0,0)exp [$ DM {R,z) -$ DM (0,0) ) + ($,(/?, z) -$„(0,0) ( p 

(3) 

where the DM potential, $ UM (R,z), satisfies the Poisson equation, 

X7 2 $ DM (R,z) = 4nGp DM (R,z). (4) 
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The solution of the coupled equations (||) and ((D for <3> DM is effected through the iterative 
scheme {n = 1, 2, 3, . . .) 

\7 2 <p n (R, z) = AnGp n ^(R, z), (5) 

where p n _i(R, z) is equal to the r.h.s. of Eq.@ with $ DM replaced by <p„^i(R, z), and 
{4>o(R, z) — <^ (0,0)} = is the initial choice for the iteration process. The quantities 
p DM (0,0) and (v 2 )*' 2 are taken as free parameters. 

Details of the iterative scheme and the numerical procedure are described elsewhere. 
After a few iterations (typically, n < 10) the potentials (p n converge towards the desired 
potential $ DM . We checked our numerical code against test equations whose exact solutions 
are known. We also check our numerical results for the actual equations and (|J) against 
analytical results for small and large values of R and z. 

Once $ DM has been calculated, the rotation curve, O(-R), is obtained through the relation 

&2{R) = (^[ $ dm(^) + M^)]) - (6) 

Note that the contribution of the visible disk to Q 2 (R) is proportional to its surface density 
[see Eq. (4-159) of Ref. 0], while that of a perfect isothermal sphere is proportional to the 
square of the velocity dispersion of its particles [see Eq.(4-127b) of Ref. ||]. 

The theoretical rotation curves thus obtained for various values of the parameters 
p DM (0,0) and (v 2 )^ 2 are to be compared with observations [|l],|Il],|l^] to ascertain the do- 



main of the parameter space which is acceptable. This comparison is shown in Fig.|] for 
S = 80 M Q pc -2 and p DM (0, 0) = 1 GeVcm -3 . The value of 80 M Q pc -2 for E , it being the 
upper limit on the allowed value of So in our calculation, gives us a conservative estimate of 
(i.e., a lower limit on) (^ 2 )™- This is because, for a given value of at a given value of R, 
a lower value of the disk surface mass density (S ) requires a higher value of {v 2 }^ 2 (for a 
fixed value of p DM (0, 0)). Our choice of p DM (0, 0) ~ 1 GeVcm -3 is dictated by the constraint 
that p DM (i? , 0) ~ 0.3 GeVcm -3 and the need to fit the rotation curve. A slightly lower 
value of p DM (0,0) generally requires higher values of (v 2 )^ 2 in order to satisfy the above 
constraint and to fit the rotation curve. In this sense, our choice of p DM (0, 0) ~ 1 GeVcm -3 
yields, again, a lower limit to (v 2 )^- A higher value of p DM (0,0), on the other hand, can 
be consistent with the constraint p DM (R ,0) ~ 0.3 GeVcm -3 for sufficiently low values of 
(f 2 ) 1 / 2 ; however j n this case, the rotation curve falls steeply beyond the solar circle and 
thus provides a poor fit to the rotation curve. 

In order to determine (a lower limit to) the best-fit value of (v 2 )^ 2 we have calcu- 
lated x 2 = E ( e ^ ) - ^ {Rl) Y as a function of (v 2 )^ 2 (for S = 80 and 4OM pc -2 and 

p DM (0,0) = 1 GeVcm -3 ), where N is the number of observational data points, Qi(Ri) and 
©i,o(-Rj) are the theoretical and observational value of the rotation speed, respectively, for 
the ith data point for which R = Ri, and Gi is the la uncertainty in the measured value 
of 0i iO (-Rj). We calculate the above \ 2 f° r the entire data set for R in the range ~ (2-20) 
kpc as well as for the restricted data set for R in the range ~ (10-20) kpc in which the 
observed rotation curve data show a conspicuous rising trend. For S = 8OM pc -2 , both 
data set give a minimum x 2 a t (v 2 ) 1 ^ 2 ~ 600 km s -1 . For So = 4OM pc -2 , the minimum of 
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the x 2 li es & t (f 2 ) 1,/2 ~ 750kms _1 for the restricted data set while the minimum is beyond 
900 km s -1 for the full data set. From the above analysis we conclude that the lower limit 
on (i; 2 ) 1 / 2 is ~ 600 kms -1 . 

\ ' DM 

Notice from Fig.[I] that for (v 2 )^ ~ 300 km s^ 1 the potential of the visible component 
concentrates the distribution of DM towards the centre, causing the rotation curve to fall 
below the observational data at large galactocentric distances. As the kinetic energy of the 
DM particles increases with increased value of (v 2 ) 1 ^ 2 the particles are affected progressively 
less by the potentials and spread out farther. This causes the rotation curves to be elevated. 

We thus see that the rms velocity of particles of DM needed to generate the observed 



rotation curve is higher than that adopted in a variety of discussions of DM |L8[ . Indeed, 



we had an inkling that this might be so, based on our analytic estimates made earlier in this 



context [19|. The implications of this result are multifarious: 



1. Since the typical velocity of individual DM particles is higher by at least a factor 
of ~ 2 on the average, the energies they would deposit in the detectors would be 
higher by at least a factor ~ 4. This would make these events stand out against the 
background. 

2. The higher velocities imply higher fluxes and the event rates would be increased by 
at least a factor of ~ 2. 

3. When the observed pulse height spectrum in the detectors are reanalysed taking the 
above two points into account the existing bounds on the masses and other properties 
of dark matter particles would become substantially more stringent. 

4. The higher velocities would also mean lower rates of capture by the Sun by accretion; 
consequently the flux of high energy neutrinos arising from their annihilations in the 
central regions of the Sun [2~U| is expected to be correspondingly smaller. 



5. The large velocities would also imply an extended halo (with an estimated mass of 
~ 1.5 x 10 12 M up to ~ 100 kpc) whose influence on the dynamical motions within our 
Galaxy and on the Local Group, as also the tidal effects on the dwarf-spheroidals would 
become important. For example, this high temperature halo will impart stability to 
the disk according to the criterion derived by Peebles and Ostriker |21[|. 



These issues are under study and will be reported elsewhere. 
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FIGURES 



FIG. 1. The theoretically calculated rotation curve of the Galaxy for various values 
°f ( v2 )dm com P are d with the available observational data [10|,ll,|l2|. All curves are for 
p DM (0, 0) = lGeVcm -3 and So = 80M Q pc -2 (see text). The data and error bars for R in 
the range ~ (2-17) kpc are from Fig. 3 of Ref. ]H| , and those for R > 17kpc are from Fig. 2 of 
Ref. £L2|. The data for R below ~ 2 kpc are from Ref. |0|. 
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1)<v 2 > l/2 DM =300 km s" 1 



2) <v 2 > 1,2 DM = 450 km s" 1 

3) <v 2 > 1/2 DM = 600 km s" 1 



4) <v 2 > 1/2 DM =750 km s" 1 

5) <v 2 > 1/2 DM = 900 km s" 1 
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